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Abstract 

We consider continuous-time average consensus dynamics in which the agents' states 
are communicated through uniform quantizers. Solutions to the resulting system are 
defined in the Krasowskii sense and are proven to converge to conditions of "practical 
consensus". To cope with undesired chattering phenomena we introduce a hysteretic 
quantizer, and we study the convergence properties of the resulting dynamics by a 
hybrid system approach. 

1 Introduction 

Communication constraints play a major role in consensus and related problems of dis- 
tributed computation and control. Such constraints can be represented by a graph of 
available communication links among agents, together with further restrictions on which 
information can be exchanged across links. Recently, the constraint of quantization, that is 
of communication restricted to a discrete set of symbols, has received significant attention. 
Although most works to-date have dealt with discrete-time dynamics, it is worth consider- 
ing the same restrictions in the context of continuous-time dynamics, because the dynamics 
of the agents is naturally described by continuous-time systems in many applications. An 
example of this is robotic networks. One might argue that it is possible to study the effect 
of quantization on continuous-time systems by considering their discretized or sampled-data 
model. However, implementing consensus control algorithms by discretizing the dynamics 
requires implicitly that all the agents sample synchronously with the same clock, a require- 
ment which is difficult to satisfy in practice. This lack of synchronicity may disrupt the 
convergence properties of the algorithm and thus asks for a different approach. This paper 
proposes an approach which deals with consensus problems in continuous-time without re- 
lying on sampled-data systems: consensus is achieved by quantized measurements which are 
transmitted asynchronously. Indeed, the quantization of the states induces a partition of 
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the space into regions and each agent transmits quantized information about its state only 
when the state crosses the boundary of the quantization regions. In this sense communica- 
tion among the agents takes place asynchronously. Throughout the paper no assumption 
is imposed on the resolution of the quantizers, which therefore can be also very coarse. 
Hence, this paper specifically aims at giving a rigorous treatment of continuous-time aver- 
age consensus dynamics with uniform quantization in communications. It is an established 
fact that consensus problems can be reformulated in terms of feedback control systems: as 
expected, when quantization enters the loop, the stabilization problem becomes more chal- 
lenging. From a mathematical point of view, a consequence of quantization is that we obtain 
a system with a discontinuous righthand side, which is not guaranteed to admit solutions in 
classical sense. This paper proves that classical or Caratheodory solutions may in fact not 
exist: considering solutions in some generalized sense is thus unavoidable. 

The literature provides different approaches to the technical problem of having systems 
with discontinuous righthand sides (see [8, 15] for a review of these topics). Tools from the 
theory of discontinuous differential equations and nonsmooth analysis have already been 
applied in consensus problems in [7, 9]. Here we focus on Krasowskii solutions essentially 
for two reasons: effectiveness and generality. With regard to effectiveness, there are many 
results available concerning the existence and continuation of Krasowskii solutions, as well 
as a complete Lyapunov theory [1, 3]. With respect to generality, since the set of Krasowskii 
solutions includes Filippov and Caratheodory solutions, results about Krasowskii solutions 
also hold for Filippov and Caratheodory solutions, provided that they exist. On the other 
hand, the set of Krasowskii solutions may be too large. In particular, from a practical point 
of view it may contain sliding modes which induce chattering phenomena. In the context 
of quantized consensus, chattering amounts to fast information transmission between the 
agents. This is undesirable because it results in algorithms which require large bandwidth 
communication channels to be implemented. To cope with this issue, we propose the use 
of a quantizer endowed with a hystcretic mechanism, and study the resulting dynamics by 
a hybrid system approach. Specifically, we provide an estimate of the data rate needed 
to implement the quantized continuous-time consensus algorithm. These results can be of 
interest to other application fields, including load balancing problems and real-time control 
systems. 

With respect to earlier literature, our contribution is twofold. On one hand, we give 
a mathematical treatment, in terms of differential equations with discontinuous righthand 
sides, of a continuous-time consensus system under uniformly quantized communication of 
the states. We do this when the communication graph is only weakly connected and weight 
balanced. After showing basic properties of solutions, such as existence, boundedness and 
average preservation, we prove convergence to a set containing the equilibria of the system. 
This set depends on the communication graph, and is reduced to the set of equilibria in the 
particular case of symmetric graphs. On the other hand, our paper is the first to propose the 
application of hysteretic quantizers to solve a consensus problem with data rate constraints. 
Preliminarly, we prove that the system is well-posed, in the sense that a solution exists, is 
forward unique, and the set of switching times is locally finite. The main results consist in 
proving convergence and estimating the required data rate. 

In our paper, we provide convergence results for both the Krasowskii and the hysteretic 
dynamics discussed above. Due to the constraint of static uniform quantization we can- 
not obtain exact consensus, but we can obtain approximations of the consensus condition 
which we informally refer to as "practical consensus" . Similar conditions have been obtained 
elsewhere in the literature, which has already considered some related problems of quan- 
tized consensus. A number of publications have discussed -mostly in discrete-time systems- 
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several options to deal with the quantization constraint, namely uniform deterministic quan- 
tizers [13, 18], uniform randomized quantizers [2], logarithmic quantizers [5], and adaptive 
quantizers [17]. Moreover, quantization has been considered also in gossip consensus al- 
gorithms [16, 6]. Regarding continuous-time systems, relevant bibliography includes [7], 
which discusses discontinuous differential equations with consensus applications, and [22], 
which studies a rendezvous algorithm in which each agent tracks another agent assigned to 
it by a quantized control law. A recent paper ([11]) is also related: using graph-theoretical 
tools, the authors study a quantized consensus problem, under the assumption that the 
communication graph is a tree. 

The present paper is organized as follows. In Section 2 we recall definitions and results 
from graph theory and about average consensus dynamics without quantization, recalling 
how the consensus problem can be reformulated in terms of a stability problem. We also 
present the problem of quantization in a continuous-time setting. In Section 3 we introduce 
state quantization and Krasowskii solutions and we study the fundamental properties of the 
system: we compute the set of equilibria, we prove average preservation, and we deduce 
asymptotic and finite-time convergence results. Then, in Section 4 we define and study 
the system under a hysteretic quantizer, in terms of a hybrid system. After proving that 
chattering can not occur, we study the fundamental properties of this system: equilibria, 
average preservation and convergence. Some simulations are given in Section 4.3 which 
illustrate our results and point to possible future research directions, which are discussed in 
our conclusions in Section 5. 

Notations. Given a subset A of the Euclidean Space M. N , we denote as A its topological 
closure, by int (A) its interior and by dA its boundary. Given N £ N, we let 1 (0) be the Ax 1 
vector whose entries are 1 (0), / be the A-dimcnsional identity matrix and f2 = I — N _1 ll* , 
where the symbol * denotes conjugate transpose. || ■ || denotes the Euclidean norm both for 
vectors and matrices. 

2 Preliminaries 

2.1 Graph theory 

Let there be a weighted (directed) graph G = (V,E,A), consisting of a node set V = 
{1, . . . , A}, an edge set E C V x V and an adjacency matrix A £ K> X N such that Aij > 
if (j, i) £ E, and Aij = if (j, i) $ E. For every node, we define the set of its in-neighbors 
as Afi = {j £ V : £ E}. We assume no self- loops in the graph, that is i £ Mi for 

every i £ V. Nodes (vertices) are referred to as agents, edges as links. Let d^ := Y^j=i Aij 
and d ou t,j = S"=i -^ij be, respectively, the in-degrec and the out-degree of node i £ V . 
A graph is said to be weight-balanced if the out-degree of each node equals its in-degree. 
Let D = diag(Al) be the diagonal matrix whose diagonal entries are the in-degrees of each 
node, which are equal to the number of incoming edges if the nonzero entries of the adjacency 
matrix A are all equal to 1. Let L = D — A be the Laplacian matrix of the graph G. Note 
that LI = 0, and that 1*L = 0* if and only if G is weight-balanced. A path in a graph is 
an ordered list of edges. Given an edge (i,j), we shall refer to i and to j as the tail and the 
head of the edge, respectively. An oriented path is an ordered list of edges such that the 
head of each edge is equal to the tail of the following one. The graph G is said to be strongly 
connected if for any i,j£V there is an oriented path from i to j in G. Instead, it is said to 
be weakly connected if for each pair of nodes i, j there exists a path which connects i and 
j. Observe that weakly connected weight-balanced graphs arc strongly connected graphs [9, 
Proposition 2]. Recall the following result, which can be derived from [4, Theorem 1.37] and 
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[9, Formula (1) and Section 2.2]. 

Lemma 1 Let G be a weighted graph and suppose it is weight-balanced and weakly connected. 
Let L be its Laplacian matrix. Then: 

(i) The matrix Sym(L) := L+ 2 L is positive semi-definite. 

(ii) Denoted by A2(Sym(L)) the smallest non-zero eigenvalue o/Sym(L), 

x*Sym(L)x > A 2 (Sym(L))||z - ~^x\\ 2 , 
for all x <G R N , where || • || denotes the Euclidean norm. 

Sometimes, it may be convenient to restrict our attention to symmetric graphs, that is 
graphs such that A = A* . For a symmetric graph there is no distinction between strong and 
weak connectedness, so in that case we shall just say that the graph is connected. 

2.2 Feedback consensus dynamics 

Let x : R>o — ► be a time dependent vector representing the agents' states. Its dynamics 
can be written in terms of the control system 

x = u, (1) 

where x, u £ M. N . Our aim is to construct a control law u : R — > K N such that, for all initial 
conditions, solutions to (1) satisfy the average consensus condition, that is 

lim x(t) = x ave (0)l, 

t— >oo 

where we let x ave (t) = N~~ 1 l*x(t). It turns out that one such control can be given in the 
feedback form u = —Lx so that the implemented system becomes 

x(t) = -Lx{t). (2) 

Componentwise, this reads as 

Xi (t) = A a ( x i (*) ~ x i (*)) » v 1 e v ( 3 ) 

where we recall Mi is the set of the neighbors of agent i. 

The following result, which can be deduced from [9], gives the weakest conditions for 
system (2) to converge to average consensus. 

Lemma 2 If the weighted graph G is weakly connected and weight-balanced, and x(t) satis- 
fies (2), then 

lim x(t) = x avc (0)l. 

t— >oc 

In view of the above result, from now on we shall make use of the following standing 
assumption, unless otherwise stated. 

Assumption 1 The communication graph G is weakly connected and weight-balanced. 
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2.3 Quantized consensus dynamics 

It is clear that the dynamics (2) is a rather idealized version of what can actually be im- 
plemented in a real control system. A very natural issue, as discussed in the introduction, 
is quantization. Let q be a quantizer, that is a map q : M — > S with S a discrete subset of 
R (If z £ M. d , q(z) is meant componentwise). Several quantized dynamics based on (2) are 
of interest from the point of view of the applications: indeed, quantization may be inherent 
either to communication or sensing among agents, or to the computation of the feedback 
control, or to its application to the system. In other words quantization can occur at the 
communication, sensor, computation or actuator level. If the sensed/communicated data 
about the neighbors states is quantized, due to the use of a digital lossless channel or to 
finite-precision sensors, the resulting dynamics is 

±{t) = -Lq(x(t)). (4) 

Despite its interest, the dynamics (4) has not received much attention yet, although its 
discrete-time counterpart has been widely studied, for instance in [18, 13]. In the present 
paper, we shall focus on this dynamics, considering two different quantizers. First, we take 
a uniform static quantizer; then, we design a hysterctic quantizer which prevents chattering. 



3 Krasowskii quantized dynamics 



Let us consider a uniform quantizer q : K - 

q(*) = 

Note that \q(z) - z\ < — . Moreover if x £ 
then consider the quantized dynamics 



AZ, defined by 
A. 



z 1 
A + 2 



pN 



, we let q(x) = (q(xi), . . . , q(xAr))*. We shall 



x = —Lq(x) 



(5) 



Since the righthand side of (5) is discontinuous, it is important to specify in which sense 
solutions have to be intended. One definition of solution is in the sense of Caratheodory : 
an absolutely continuous function x(t) is a Caratheodory solution to (5) on an interval I if 
it satisfies (5) for almost every t £ I. Such solution is said to be complete if / = [0, +oo). 
However, this natural definition is not suitable to study the system at hand, because of the 
following fact: 



Proposition 1 (Caratheodory solutions) There are weakly connected and weight-balanced 
graphs and initial conditions for which no Caratheodory solution to (5) exists. Moreover, 
in these cases, there is a positive-measure set of initial conditions such that Caratheodory 
solutions starting from that set are not complete. 

Proof: The righthand side of (5) is constant over open hypcrcubes with edge of length 
A, and surfaces of discontinuity are hyperplanes. Each hypcrplane is orthogonal to one of 
the elements of the canonical basis of R , ex, ■ ■ ■ , &n- Let us consider one such hypcrplane, 
namely one which is orthogonal to ej and let us denote it by Sj. Let us consider a point 
x £ Sj, i.e. such that Xj = (k + i)A for some k £ Z, and such that if i ^ j, then 
%i {h + -|)A for any feeZ. Let I(x) be a neighborhood of x such that, for any x £ I(x), 
it holds that Xj ^ (h + i) A for any integer h ^ k, and Xi ^ (h+ i) A if i ^ j for any h £ Z. 
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Let /+(£) = {x E I(x) : Xj > (k + i)A} and I~(x) = {x E I(x) : xj < (k + ±)A}. Let q 
be the value of q(x) in I + (x), f + be the value of —Lq(x) in and /~ be the value of 

— Lq(x) in I~(x). Then, letting the nonzero entries of the adjacency matrix to be all equal 
to 1, the components / + and / _ which are orthogonal to Sj have simple expressions and 
interpretation: 

ft : = (f + , e j) = fj~ = ^2 9k- dinjqj, 
/J:= (f-,e j ) = fr = f+ + d in , j A, 

where d m j is the in-degree of node j. If either ft is positive or /J < — d- m jA, then 
and /JJ have the same sign, otherwise they have opposite sign. A numerical example for 
the latter case is as follows: let d- m< j = 2, qj = 2A and {qk ■ k E Afj} = {A,2A}, then 
= —A and fj = A. In this case, /J < and /J > 0, so that there are no Caratheodory 
solutions starting at x. Moreover, solutions starting in a (sufficiently small) neighborhood 
of x reach the surface in finite time, and can not be extended further. ■ 

Example 1 Consider a graph with adjacency matrix whose entries are: A\2 = ^21 = 
1,^23 = A32 = 1,^4.11 = A22 = ^33 = A13 = A31 = and the initial condition x — 
(A, |A, 2A). Clearly the graph satisfies Assumption 1. Note that the righthand side of (5) 
is discontinuous at x, being on the plane X2 = § A. Following the proof of Proposit ion 1 one 
can easily check that f% = —A and f% = A. A possible Caratheodory solution issuing from 
x could not leave the surface X2 = §A in the direction of decreasing X2 because f% — A > 
nor could leave the surface in the direction of increasing x%, since /+ = -A < 0. The 
only possibility for a solution would be to remain on the surface but such a solution would 
not satisfy the definition of a Caratheodory solution. Hence, there are no Caratheodory 
solutions issuing from x. 

In view of Proposition 1, in the sequel of this paper we shall consider solutions in a 
more general sense, which is due to Krasowskii. An absolutely continuous function x(t) is a 
Krasowskii solution to (5) on an interval / if it satisfies at almost every t £ I the differential 
inclusion 

x G IC(-Lq(x)), (6) 

where 

IC(-Lq(x))= f]co(-Lq(B(xJ))), 

8>0 

and B(x,6) is the Euclidean ball of radius S centered in x. The solution is said to be 
complete if I = [0, +00). Note that, thanks to Theorem 1 in [20], we have 

K{-Lq(x)) = -LK{q{x)) C -L( X /C(q(x,)) = X ( V A tJ (/C(q(^)) - JC(q( Xl )))) (7) 

where X j denotes the Cartesian product of the sets indexed by i. Representations of the 
map q and of the set-valued map JCq arc given in Figure 1. 

Lemma 3 (Krasowskii solutions) For any x° € ~BL N there exists a complete Krasowskii 
solution x(t) to (5) such that x(0) = x . 



G 



Figure 1: Visualization of the map q(x) and the set- valued map Kq(x), when A = 1. 



Proof: Local existence of Krasowskii solutions is guaranteed by the fact that the 
righthand side of (5) is measurable and locally bounded (see [15]). By standard arguments 
completeness of solutions can be deduced by their boundedness (see, e.g., [19]). We then 
prove that solutions are bounded. Let x(t) be a Krasowskii solution to (5) such that x(0) = 
x°. Let m(t) = mm{xi(t),i e V} and M(t) = max{xi(i),i £ V}. Let t > be fixed. We 
prove that m(t) > m(to) for all t >t . Analogously it can be proven that M(t) < M(to) for 
all t > to- By taking to = we will then deduce that for any t > and for any i £ V one 
has m(0) < x%{t) < M(0), i.e. x(t) is bounded. We first consider the case m(to) (k ~ |)A 
for any fc £ Z. Assume by contradiction that there exists i > to such that m(t) < m(to) 
and let i € V be such that xj(t) = m(t). Note that m(t) is a continuous function since it 
is the minimum of a finite number of continuous functions. Then there exists S > such 
that m(t — 5) = m(t ) and for all t £ (i — 6,t) one has q(m(t)) = q(xj(t)) = q(m(io)) 
and q(xi(t)) > q(m(to)) for all i E V. We remark that for all t £ (i — S,t) it holds 
K.(q(xj(t))) = {q(m(to))} and v > q(m(t )) for all v £ K,(q(xj(t))) and any j £ V. Since 
±i(t) £ (-L^(q(x(t)))i C EiGWi^iW^W)-^^*))). we § et that ^ for 

almost all t £ (t — 5,t). On the other hand, xj(t) = m(t) < m(to) = m(t — 5) < xj(t — 5). 
Then there must exist a subset / of (i — S, t) such that / has positive measure and ij(t) < 
for almost all t £ I, i.e. a contradiction. We now consider the case in which m(t ) = (k— |)A 
for some k £ 7L. Let t* = mi{t > to : m(t) < m(to)}, and assume by contradiction that 
t* < +oo. Note that m(t) > m(to) for t £ [to,t*) and, since m is continuous, one also 
has m(t*) = m(to). Let us fix any t > t* such that (k — |)A < m(t) < m(to) for any 
t £ (t*,t). We take a sequence {£n}neN such that t n £ (t*,t) for any n and t n —> t* . Since 
m(t n ) < m(to), thanks to the previously analyzed case, we get that m(t) > m(t n ) for all 
t > t n . By passing to the limit we get that m(t) > m(t*) = m(to) for all t > t* and finally 
we have that m(t) > m(to) for all t > to, again a contradiction. B 

We observe that Lemma 3 does not guarantee uniqueness of solutions, as explained in 
the following remark. 

Remark 1 (Sliding mode) From the proof of Proposition 1, we see that if Carathcodory 
solutions originate from any point x £ Sj such that f± < and /J > 0, then they do 
not exist, while Krasowkii solutions do. Krasowskii solutions starting in a neighborhood of 
x reach the surface in finite time and slide on it. We remark that such solutions are not 
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backward unique. A general discussion about these behaviors can be found in [12, page 51]. 

Example 1 (Cont'd) Consider again the system in Example 1. It is straightforward to 
check that 







1 ■ 


i 


o / 







)C(-Lq(x)) = col I -A 1,1 A I U^el 3 :t.= | (1 - 2A)A I , A e [0, 1 
IV °/ V- A /J I V-(1-A)A/ 

Hence a Krasowskii solution issuing from x is the solution to x\ = -j , x% = 0, ±3 = — % 
with initial condition x. On the other hand, as explained previously in Example 1, no 
Caratheodory solution issuing from x exists. 

A preliminary result establishes that the average of the states is preserved by Krasowskii 
solutions to (5): this fact will be a key step to obtain many of the following results. We 
recall that this result holds under the standing Assumption 1: the same is true for all the 
following ones, with the exception of Theorem 2 and Proposition 4 which require a stronger 
assumption. 

Lemma 4 (Average preservation) Letx(t) be a Krasowskii solution to (5). Then, x ave (i) 
a^ave(O) for all t > 0. 

Proof: By definition, x avc (t) = N~ 1 l*x{t), and then ± ave (i) = N^^ifr). By (6), 
iavoW E iV _1 l*X:(-iq(x(*))) fora.e.t. 
Thanks to the first equality in (7) and the fact that 1*L = 0*, 

iavcW E -N-H*LIC( q (x(t))) = {0}, 
then we get that x ave (t) = for a.c. t and finally that £ ave (f) is constant. ■ 



3.1 Graph-dependent convergence results 

A first set of results regards the limit behavior of system (5), which depends on the quantizer 
and the graph topology. The following proposition proves convergence of solutions to a 
certain set. 

Theorem 1 (Convergence) If x(t) is any Krasowskii solution to (5) and 
M^R»:^||*-W0)1|I< ^|_f}, 

then dist(ai(i), M) ^0 as t +00. 

Proof: Let y(t) = Qx(t) = x(t) - x ave (t)l. Then, y = fix e MC(-Lq(x)). Since 
/C(Q(— Lq[x))) = fllC(— Lq[x)) = K,(—Lq(x)) = —LK.q(x), we have y € — LJCq(x). Consider 
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the function V(y) = \y*y and let v G K,q(x). Note that if v G 1Cq(x), then ||« — x|| < V^Vf- 
We have that 

W(y) • y = - y*Lv 

= — y*L(x + v — x) 
= — y*Lx — y*L(v — x) 
= - y*Ly-y*L(v - x) 
= - y*Sym(L)y - y*L(v - x) 

<-A 2 (Sym(L))|M| 2 + ||y||||L||^ViV 

= -A a( g^(L))||y||(||y||- ] ^^^ 

where the inequality in the second-last line follows from Lemma 1. This implies convergence 
to the set 

Finally the statement follows from average preservation in Lemma 4. ■ 

With a slight extension of the argument leading to Theorem 1, we can prove finite-time 
convergence to a set larger than the set M in Theorem 1, and provide an estimate of the 
convergence speed. 

Corollary 1 (Finite-time convergence) If x(t) is any Krasowskii solution to (5), then 
for any e G (0, 1) there exists a finite time T(e) such that x(t) belongs to the set 

for all t > T(e). 

Proof: Consider, as in the proof of Theorem 1, the differential inclusion y G —LKq(x), 
the function V(y) = \y*y and let v G /Cq(x). Then, as before, 

VV(y) -y<- A 2 (Sy m (L))||2/|| (\\y\\ - —B—^Vn 



A 2 (Sym(L)) 2 



From the latter, we see that, if 



then 



VV(y)-y<-e\ 2 (Sym(L))\\y\\ 2 
= -2e\ 2 (Sym(L))V(y) 

If at some time to the condition (8) is satisfied by y(to), then there exists t > to such 
that the Lyapunov function computed along the trajectories of (5) satisfies 

V(y(t)) < e - 2eA =( s y in ( i )))(*-*o)^(y(t )) 
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and therefore 

\\y(t)\\<e- eX2iSym{L)mt - to) \\y(to)\\- 

Assuming without loss of generality that ||y(io)|| 7^ 0, from the latter inequality we conclude 
that, if at time to the condition (8) is satisfied, then there exists a time 

f "I ( 1 llXrll A y/~N \ 1 

T{£) = ^ °' eHSym(L)) Hni»)jTil (9) 



such that ||y(i)|| satisfies 



for any t > to + T(e). The thesis then follows, recalling the definition of y(t) and Lemma 4. 

■ 

The above convergence results assert that the error induced by quantization, with respect 
to the non-quantized consensus dynamics (2), can be made arbitrarily small by decreasing 
the quantization error. Up to a \f~N factor, due to the length of the vector, the committed 



error is proportional to and then depends on the network. It is plain that it would 

be of interest to state a result of convergence to a stronger practical consensus condition, 
in which the committed error does not depend on the network topology but only on the 
quantizer precision. This issue is the topic of the next paragraph. 



3.2 Equilibria 

In this paragraph, we shall describe the equilibria of the system (5), which depend on the 
quantizer precision only. Hence, proving convergence to equilibria turns out to be a way to 
prove a practical consensus condition which does not depend on the network. The following 
proposition characterizes the equilibria of the system. We recall that x is a (Krasowskii) 
equilibrium if the function x(t) = Xq is a (Krasowskii) solution, that is if G /C(— Lq(xo))- 
Let 

V = {x G R N : 3k G Z such that q(xi) = Afc, Vi G V}. 

Proposition 2 (Equilibria) The set of Krasowskii equilibria of (5) is T>. 

Proof: Note that Lq(x) is zero in T>, and is discontinuous on the boundary of T>. 
Let us define the set of Krasowskii equilibria as E — {x G R N : G — LJC(q(x))}, and let 
E = {x G R N : G —L(XiJC(q(xi)))}, where we recall that Xj denotes the Cartesian product 
of the sets indexed by i. Since Theorem 1 in [20] implies that JC(q(x)) C Xj/C(q(xi)), wc 
have that E C E. In order to prove that T> = E, we will first prove that int (2?) C E. Since 
E is closed due to the fact that the set-valued map K{q{-)) is upper semicontinuous (see, 
e.g., the definition of upper semi- continuity given in [10]), then also T> C E C E. Later we 
will prove that ECU. These two facts imply that T> = E = E, and namely our statement. 
Let us prove that T> C E. Let us assume xo G int 2?. Since q is continuous at xq, then 
fC(q(x )) = {q(x )} and K(-Lq(x )) = {-Lq(x )} = {-L(Akl)} = {-AkLl} = {0}, i.e. 
xq G E. The points x G V also belong to E thanks to the fact that E is closed. Then, 
let us prove that E CP. xo G E if there exists v G X,/C(q(a;o i )) such that Lv = 0. This 
is equivalent to the fact that kerL n x i IC(q(xo i )) ^ 0. Since kcrL = spanl, there exists 
V G kcrL n x i K,(q(xo i )) if there exists A £ 1 such that v = XI, i.e. Vi — A for any i G V 
and Vi G K.(q(xo i )). Such A can be either A = Afc for some fc G Z, or A ^ Afc for any k G Z. 
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In the first case, for any i £ V we have that {xi : A G /C(q(a;,))} = \(k — ^)A, (k + |)A] . 
In the second case, we have that {xi : A G /C(q(xi))} = {([-^J + l)A}- Finally we get that 
if v = Al € x j/C(q(a;o i )) then for every i € V either £o; G \{k — ^)A, (fc + ^)A] for some 
£ Z, or = (k + ^)A for some k G Z, i.e. io £ D. I 

We now prove that T> is strongly invariant, i.e. there are no trajectories exiting T>. 

Proposition 3 (Strong invariance) // x(t) is a Krasowskii solution to (5) such that 
x(0) € V, then x(t) G V for all t > 0. 

Proof: Recall that 

V = {x G R N : 3k G Z s. t. (fc - i)A < a;, < (k + i)A, Vi G V}. 

Let x(i) be a solution to (5) such that x(0) = x° G V. If x° G int(X>), then Lq(x°) = 0, 
and thus K,Lq[x°) = {0}. Hence int (T>) is invariant. Let then x° G 9P. Then there 
exist fc° G Z and C V (not both empty) such that x;(0) = (fc - 1/2)A for all 

i G V - , arj(O) = (fc° + 1/2)A for all i G V+, and ^(0) G ((fc° - 1/2)A, (fc° + 1/2)A) 
for all i € V \ (V - U V + ). Let us assume by contradiction that there exists T > such 
that ie(T) ^ X>, i.e. either there exists ! £ V such that xj(T) > (k° + 1/2) A or there 
exists j G V such that xj(T) < (k° — 1/2) A. For brevity, we examine only the former 
case. Let i* be such that Xi*(T) = max{.T.;(T), i G V}. Without loss of generality, we can 
assume that (k° + 3/2)A > Xi *(T) > (k° + 1/2) A and x,(T) < (k° + 3/2)A for all i G V. 
Thanks to the continuity of x(t), there exists T 1 < T such that x 4 »(T') = (k° + 1/2) A and 
(fc° + 1/2)A < ^.(i) < (fc° + 3/2)A for all t G (T',T). Since x^(T) - ^.(T') > there 
exists a subset T of (T',T) such that T has positive Lebesgue measure and for all f G T 
the derivative Xi*(t) exists and is positive. On the other hand for all t G (T',T) one has 
(fc° + 1/2) A < (t) < (fc° + 3/2)A, which implies q(x;. (i)) = (fc° + 1)A for all t G (T', T). 
Let us now consider any v G K,(q(x(t)) with £ G (T',T). It holds Uj» = (fc° + 1)A and 
Vj < (k° + 1)A for all j^i*, and therefore (-Lu) 4 . < d in ^ (fc° + 1)A - d hl ^ (fc° + 1)A = 0. 
From this fact it follows that Xi*(t) < 0, i.e. a contradiction. ■ 

Next, we provide a second convergence result, stating that on any connected symmetric 
graph the quantized dynamics converges to the set of equilibria. This fact implies that the 
error induced by quantization does not depend on the network properties, but only on the 
quantizer. 

Theorem 2 (Convergence to equilibria) If the weighted graph G is symmetric and x(t) 
is any Krasowskii solution to (5) ; then dist(x(t), T>) — > as t — > oo. 

Proof: We consider the function V(x) = ^x*x and we prove that 

VV(x) • (-Lv) < 
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for every x and for every v £ ICq(x). Let us first remark that ICq(x) C Xj/C(q(xj) and then 
v £ JCq(x) implies that i> £ X j/C(q(:Ei)). 

W(x) • (--Lw) = - x*Lv 

i 3 

= ~ 2 ^2( X 3 ~ X i) A ij( v 3 -Vi)<0, 
ij 

where in the third equality we have used the symmetry of the graph, i.e. of the matrix A, 
and where the last inequality is a consequence of the fact that if a, b £ M, a < b, then for 
all u £ /Cq(a) and all w £ ICq(b) it holds u < w. Namely, all terms of the above summation 
are nonnegativc. Let Z = {x £ l w : 3v £ JCq(x) such that W(x) ■ (—Lv) = 0} and let M 
be the largest weakly invariant subset of Z. By LaSalle invariance principle for differential 
inclusions (see [3], Theorem 3) we have that any solution x{i) to x £ — L(JC(q(x)) is such 
that dist(x(t), M) as t -> +oo. 

We now prove that ZCP. More precisely we prove that if x ^ T> then x ^ Z. Proving 
x g" Z is equivalent to prove that for all v £ /Cq(S) one has W(x) ■ (—Lv) < 0, i.e. 
Sijfe — %i)Aij(vj — Uj) > 0, and, for this to hold, it is sufficient to prove that there exist 
i, j such that A^ ^ and (xj — Xj){vj — Vi) 7^ 0. Let us then assume that x ^ T>. 

This means that for all k £ Z there exists i G V such that xj < (k — 1/2) A or there exists 
j £ V such that x j > (k + 1/2) A. Let k\ = q(^i). Then there exists i £ V (without loss 
generality we assume i = 2) such that X2 < (ki — 1/2) A or X2 > (fci + 1/2) A. We consider 
the case X2 < (fci — 1/2) A. Let ki be such that q(x2) = ^A. Clearly k2 < fci. We examine 
the following possible cases: 

(a) xx ^ (fci - 1/2)A and x 2 + (k 2 - 1/2) A; 

(b) xi = (fci - 1/2) A and x 2 ^ (k 2 - 1/2) A; 

(c) x x ^ (fci - 1/2) A and x 2 = (k 2 - 1/2) A; 

(d) xx = (fci - 1/2) A and x 2 = (k 2 - 1/2) A. 

In case (a) for all v £ /Cq(x) one has v\ = fci A and V2 = fc2A then v% — V2 = A(k% — k 2 ) 7^ 
since fci k 2 . In case (b) we remark that, since x g" T>, then k 2 < ki — 2. For all i> £ /Cq(x), 
one has ui = afciA + (1 — a)(ki — 1)A = A(fci — 1 + a) with a £ [0, 1], and v 2 = k 2 A. 
Then v\ — v 2 = A (fci — fc'2 — 1 + a) > A ^ 0. Analogously in case (c) we get v\ = Afci, 
v 2 = a(k 2 — 1)A + (1 — a)k 2 A, then V\ — v 2 = A (fci — fc 2 + a) with a £ [0,1], then v± — v 2 7^ 0. 
Finally in case (d) we also have fc2 < fci — 2, and moreover, since ICq(x) C X,/C(q(aij)), for 
all v £ K.q(x) we have 

ui =A[ai(fci - 1) + (1 - ai)fci] 
v 2 =A[a 2 (k 2 - 1) + (1 - a 2 )k 2 ], 

with ai £ [0, 1], i = 1, 2. Then we get ui — U2 = A[ki — k 2 — ol\ + CI2] 7^ 0. Hence we have 
proved that Z C D. This fact also implies that Z C D, and finally, dist(x(t),2?) — > as 
t -> +00. ■ 

We remark that Theorem 2 does not imply that solutions converge to a point in T>. 
However, Lemma 4 implies that solutions whose initial conditions belong to the hypcrplanc 
J2iLi x i — N(k + |)A for some k £ Z converge to the point (fc + ^)Al. 
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One may ask whether solutions to (5) reach T> in finite time. This claim is supported by 
numerical simulations, which we will illustrate in Section 4.3. Moreover, the following result 
shows that the claim is true for almost any initial condition. 

Proposition 4 (Finite-time convergence to equilibria) If the weighted graph G is sym- 
metric and x(t) is any Krasowskii solution to (5) such that x aV c(0) ^ {k + 1/2) A for every 
fceZ, then there exists T such that x(T) £ T>. 

Proof: Let us consider any Krasowskii solution x(t) to (5) such that a;(0) T> and its 
composition with the function V (x) = \x*x. For almost all t > we have that x(t) = —Lv(t) 
with v(t) £ ICq(x(t)) and (see the proof of Theorem 2) as long as x{t) £" T>, we also have 
that for every t > 0, 

jV(x(t))=VV(x(t))-x(t) 

= -\ X>i(*) - *i(*))>M«i(*) - 

ij 

<-\ay< fa (*)-**(*))(«*(*) -**(*)) 

(i,j)GB 

= -^X lafj (*) — ac< (*) 1 1 "rf (*) — «*(*)! 

(i,j)GB 

< - , P 1 ^ 1^' (*) _ X i (*) i I U J (*) _ U * (*) I 
1 

(i,7)TB ' 



< — - AA max |xj (t) — Xi(t)\ 



< - 0AA, 



1 

I'.Jev 1 " 1 '-' ' '~ lJ ' T 2 t>b L (i,7)Ts 1 

> 0, by integrating over the interval [0,T], we get 



where A = ruin {Ajj : Ay ^ 0} and </> = - inf { max \ Xj (t) -Xi(t)\ : x(t) € R N \T>}. If 



V(x(r)) - V(x(0)) < -c/)AAT. 

Assume by contradiction that x(t) does not reach T> in finite time. Then by letting T —¥ +oo 
we get that limT-s.+oo V(x(T)) = — oo, which contradicts the fact that V is positive definite. 

We now prove that if <f> = 0, then a; avo (0) = (k + 1/2)A for some fc £ Z. If = 0, then 
for any n € N there exists f„ such that x(t n ) £ M. N \ V and, for all £ E, it holds 

\xj(t n ) — Xi(t n )\ < — . Let us fix any ieV, and consider a path of length K connecting all 
the vertices of G and starting from k. Let us denote the states of the nodes in this path as 
X\, ...,xk, noting that some of the nodes may appear in the list more than once. From the 
fact that cf> = we deduce that 

Xl{t n ) - - < X 2 {t n ) < Xi{t n ) + -, 

n n 

2 1 12 

xi(t n ) < x 2 {t n ) < x 3 (t n ) < x 2 (t n ) + - < x 1 (t n ) + -, 

n n n n 

etc. 

- , \ K ~ , \ ~ f s K 

X\{t n ) < X K {t n ) < Xi(t n ) H . 

n n 
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By summing N of these inequalities corresponding to the N different nodes we get, since by 
definition x\ = Xk, that 



K 



N 



K 



Xk y<r t 



i—l i—1 i—1 

Since ± Y,iLi %i(t n ) = jf Y,iLi we obtain 

1 N 

— ^x 4 (0) - x k (t 



< 



K(K + 1) 



2n 



and then Xk(t n ) — > Yli=i x i(9) as n ~ * +°°- By the arbitrariness of k, we get x(t„) — > 
x avo (0)l asn-> +00. We recall that x(t n ) G M. N \ V for every n, then x aV c(0)l G R N \ V. 
Finally since \V n spanl = {(fc + |)A : fc G Z} we get that a; avc (0) = (k + §)A for 
some fc G Z. ■ 



Remark 2 (Comparison with discrete-time consensus) The results of "practical con- 
sensus" in Theorem 1 and 2 can be related with results about discrete-time consensus sys- 
tems, as those in [13] and [16]. In [13], the authors consider x(t + l) = x(t) — ^ — ^—^Lq(x(t)), 
where d ma x is the largest in-degree in G, and e > 0. Clearly V is the set of equilibria for 
this system, but its limit behavior shows limit cycles which are not contained in the closure 
of T>. However, for some example topologies (rings and complete graphs) the system can 
be proved to approach average consensus up to the quantizer precision. The paper [16] 
considers a discrete-state dynamics in which the agents communicate in randomly chosen 
pairs, and proves convergence to a set of the form {x G Z N : i, g {L,L+ € Z}. 



4 Chattering-free quantized dynamics 

The analysis in Remark 1 has pointed out the existence of sliding modes in the system 
x = — Lq(x), a phenomenon which is not acceptable in practical implementation. In this 
section we discuss a different quantization scheme to overcome this difficulty and we analyze 
the resulting system. 



4.1 Hysteretic quantizer: hybrid model 

The different quantization scheme is based on a quantizer with hysteresis which we introduce 
below. Let us consider the multi-valued map qh defined as (Figure 2) 

qhlrj ~\ jA+f jA<r<A + jA, jgZ [ } 

for the scalar r, with the understanding that if x G K™, q_h(x) = (qh(a^i), ■ • ■ , qh(a^n))* The 
evolution of qh(x(t)) as a function of x(t) can be described as follows. At time t = 0, 
qh(x(0)) = q(.x(0)). Let, for the sake of notational simplicity, qh := q^(x(t)), and qj|" := 
lim s ^ t+ qh(z(s)), x := x(t). If q h G (AZ U (AZ + A/2)) N , then 

q h + # ifa;>qh + f 

' qh-f ifx<q h -# (11) 

q h if qh - — < x < qh + -■ 
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Figure 2: The uniform quantizer with hysteresis qh(z), when A = 1. 



The rule (11) is illustrated in Figure 3. 

Suppose that in the dynamics (4), each agent quantizes the information using q^ rather 
than q. This leads to a system which can be better described and analyzed using the 
formalism of hybrid systems. To this end, we adopt the notations in [14]. Let q <G (AZ U 
(AZ + A/2)) N be the discrete state, x G R N the continuous state and X = R N x (AZ U 
(AZ + A/2)) N the state space where the system evolves. The continuous dynamics of the 
system are described by 

• : o i5 < i2 > 

which are valid as far as the state (x, q) belongs to the subset of the state space: 

A A 

C = {(x,q) G X : Vi G V,-— + q t < x, t < — + 

If on the other hand (x, q) belongs to the set 

D = {(x, q) G X : 3 i G V, such that xi < —— + qi or xi > — + 
then the following discrete update occurs: 



X + = X 












A 

2 


if Xi > — + qi 




it = < 




A 

2 


iixi < -— + qi 


(13) 




k Qi 




otherwise. 





Observe that CUfl = I and C n D = 0. The equations (12) and (13) together with the 
sets C, D define the hybrid model associated with the multi-agent system in the presence of 
the quantizers (10). In what follows we let z = (x* q*)* be the entire state of the hybrid 
model and f(z), g(z) the maps on the righthand side of (12) and, respectively, (13). Hence, 
the system is concisely described as 



i = /(*) zeC 

= g(z) zeD 



(14) 
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Figure 3: The figure illustrates the evolution of q^(x) as a function of x as described in 
(11). Suppose initially that qh — y < x < qh + y and that qh(x) = qh. As x evolves, the 
point (x, qh(x)) lies on the segment CD, and it may hit the point C or the point D, thus 
triggering a discrete transition. If the former occurs, the discrete value takes a new value, 
namely qj|~ = qh — y, the point (x, qh(x)) jumps to G and can then move towards E or 
in the opposite direction. On the other hand, if (x, qh(x)) hits D, then the quantization 
level takes the new value qh + y, i.e. qj^ = qh + y- ^ n ^ ne S ra P n ' this corresponds to the 
transition from point D to point F. The state x can then further increase or decrease. 
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with an initial condition belonging to 

X Q = {(x, q) £ X : q t - y < x t < q t + y , Vi £ V}. 

Note that initial conditions of the form (xo,q(cco)) , with xq £ R , belong to this set. 

We recall from [14] the notion of hybrid time domain and solution for a hybrid system. A 
hybrid time domain is a subset of R>o x N which is the union of infinitely many intervals of 
the form [tj, tj + {\ x {j}, where = to < t\ < ti < . . ., or of finitely many such intervals with 
the last one possibly of the form [tj, ij+i] x {j}, \tj,tj + \) x {j}, or [tj, +00) x {j}. Let z(t,j) 
be a function defined on a hybrid time domain dom2 such that for each fixed j, 1 1— > z(t,j) 
is a locally absolutely continuous function on the interval Ij = {t : (t,j) £ domz}. The 
function z(t,j) is a solution to the hybrid system (14) if z(0, 0) £ C U D and the following 
conditions arc satisfied: 

• For each j such that Ij has non-empty interior, 

z(t,j) = f(z(t,j)) fora.e.ieij 

z(t,j) £ C for all i G [min ij , sup ij ) 

• For each (t,j) £ domz such that (t, j + 1) G domz, 

<i,j + l) = j)) 
«(t,i) 6 ^. 

The solution z is nontrivial if domz contains at least another point different from (0, 0) and 
is complete if domz is unbounded. 

4.2 Hybrid model analysis 

In this subsection we detail the analysis of system (14). After proving basic properties about 
existence, uniqueness and completeness of solutions, we verify that chattering can not occur. 
We also compute the equilibria of the system, and present a convergence result. 

Lemma 5 (Basic properties of solutions) For each z(0, 0) € X$, there exists a non- 
trivial solution z(t,j) to (14), cind such solution is forward unique and complete. Moreover, 
every interval Ij = {t : (t,j) G domz}, possibly with the exception of Iq, has a non-void 
interior. 

Proof: There are two cases: z(0, 0) e D or 2(0, 0) € C. If z(0, 0) g D, then 2(0, 1) € C. 
Indeed, let i € V be any index such that Xi(0, 0) = y +<Zi(0, 0). Then Xi(0, 1) = Xi(0, 0) and 
<?i(0, 1) = y + qi(0,0), from which Xi(0, 1) = qi(0, 1). Similarly, for any index k for which 
x k (0,0) = + q k (0, 0), we have x fc (0, 1) = x fe (0, 0) and q k (0, 1) = -f + ^(0,0). For all 
the remaining indices ^(0,0) - y < xt(0, 0) < <^(0, 0) + f- and then <^(0, 1) = qe(0,0) 
and x^(0, 1) = xg(0,0). Then z(0,l) € C and the discrete transition stops in one step. 
Since the righthand side of (12) is constant, the solution z(t,j) starting either from z(0, 1) 
(if 2(0,0) € D) or from z(0, 0) (if 2(0,0) € C) exists for t sufficiently close to and is 
unique. As a matter of fact, one can compute explicitly the solution z(t,j) on the interval 
Ij, with j = (if 2(0, 0) £ C) or j = 1 (if 2(0, 0) £ D). Define the auxiliary Cauchy problem 
C(t) = —Lq(0,j), C(0) = x(0, j), where by a slight abuse of notation we are denoting the 
time variable in the auxiliary system by the same symbol t which appears in the hybrid 
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time (t,j). The problem admits the unique solution which is defined for all t > 0. Let 
t be the minimal time greater than zero for which there exists an index k G V such that 
= — + qk(0,j) or Cfc(*) = y + Qk{0,j)- Set = [0,i] if a finite t exists, otherwise 
set Ij = [0, +00). Finally let z(t,j) — (£(t)* q(0,j)*)* for all t G Ij. Then, by construction, 
z(t,j) satisfies the condition (12). If sup/j = +00, we have concluded the argument. If 
instead sup/., < +00, z(t,j) satisfies the differential equation z(t,j) — f(z{t,j)) = —Lq(0,j) 
for all t G Ij except at t = sup/j; but z(swplj,j) G D and at time t = supij a discrete 
transition occurs. After this transition, the solution can be uniquely extended starting from 
z(suplj,j + 1) G C, repeating the argument which have been just discussed. Then, z(t,j) 
is forward unique because it results from the concatenation of unique solutions to (12) and 
(13). 

Let us show that the solution is complete. By contradiction, domz is bounded, that 
is domz is the union of finitely many intervals of the form [tj,tj + {\ x {j}, with the last 
interval either of the form [tj,tj + i] x {j} or [tj, ij+i) x {j} and tj + i < +00. This yields a 
contradiction, because in view of the form of the last interval, no more discrete transitions 
take place and by (12) the continuous evolution exists for all t. 

Finally, the interval Ij = {t : (t,j) G domz}, with j such that z(0,j) G C, has a non-void 
interior because the solution to (14) from z(0,j) flows continuously with bounded velocity. 
If supij = +00, then (Ij,j) is the last interval of dom z and the proof ends. Otherwise, 
if sup Ij < +00, then by definition z(supij,j) G D and there is an index i G V such that 
either x%(suplj,j) = y + <7i(supij, j) or Xi(sup Ij, j) = — 4 + qi(suplj, j). Then, as in the 
first part of the proof, one proves that z(swplj,j + 1) G C and again Ij + i has a non-empty 
interior because the system flows continuously with bounded velocity. The thesis descends 
by induction. ■ 

The following lemma proves boundedness of solutions, a property which is useful in order 
to prove that chattering does not occur. 

Lemma 6 (Bounded solutions) For every initial condition z(0, 0) G Xq, the solution 
z(t, j) satisfies 

-f + min{ gi (0,0), . . .,q N (0, 0)} < ?i (t,j) < max{ 9l (0, 0), . . .,q N (0, 0)} + f 
-A + min{ gi (0, 0), . . . , q N (0, 0)} < Xi (t,j) < max{ gi (0, 0), . . . , q N (0, 0)} + A , 

for all i G V and for all (t,j) G domz. 

Proof: Consider the function U(z) = Taax{q±, . . . ,q^} and compute it along the 
solution to (14). During continuous evolutions, i.e. if z(t,j) G C, U(z(t)) remains constant 
since, by the second equation in (12), no component of q(t,j) changes its value. If z(t,j) G D 
and (t,j) (0,0), then U(z(t,j + 1)) < U(z(t,j)). In fact, suppose this is not true. Then, 
since q(t,j) G (AZU (AZ + A/2)) N , and at the transition each component qi either remains 
the same or changes its value by ±A/2, i.e. qi(t,j + l) = qi(t,j) or qt(t,j+l) = qi(t, j)± A/2, 
the only possibility for having U(z(t,j + 1)) > U(z(t,j)) is that one of the agents, say 
k, for which its discrete state qk(t,j) equals max{qi(t,j), . . . , qn{t, j)}, increases its value 
of +A/2 . But for this to occur, it must be true that during the continuous evolution 
which preceded the transition, the continuous state Xk(t,j) has increased its value until 
it reached the threshold y + qk(t,j) = Xk(t,j). This is a contradiction because, since 

qk(t,j) = max{<?i(£, j), . . .,q N {t,j)}, in the equation x k (t,j) = J2f=i ^y(?i(*>i) - Qk(t,j)) 
the righthand side is non-positive. Hence, we conclude that U{z(t,j)) < U(z(0,k)) for 
all (t,j) G domz such that t > and j > k, where k G {0,1} is the smallest integer 
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such that z(0,fc) £ C. In case k = 0, then U(z(t,j)) < U(z(0, 0)). In case fc = 1, then 
t/(z(0,l)) < U(z(0,0)) + A/2, and U(z(t,j)) < <7(z(0, 1)) < (7(z(0,0)) + A/2. Similarly 
one can prove that J^(z) = min{gi, . . . ,qN} satisfies W(z(t,j)) > 1^(^(0,0)) — A/2. 
The conclusion on the continuous state follows easily since z(0, 0) £ Xo implies that for each 
(t, j) £ domz, for each i £ V, qi(t,j) - f < %i(t,j) < qi{t,j) + f-. ■ 

The last part of the thesis of Lemma 5 implies that the set {t £ M.>o\z(t, j) £ D} has a 
zero measure. Moreover, a finer analysis can be performed. Let us consider the switching 
times of a solution z(t,j) to the hybrid system originating from z(0, 0) £ Xo, i.e. the times 
tj such that (tj,j) £ domz implies (tj,j + 1) £ domz. The set of switching times is locally 
finite if, for any compact subset of M>o, there is only a finite number of switching times 
which belong to that compact subset. 

Proposition 5 (No chattering) For every initial condition z(0, 0) £ Xq, the set of switch- 
ing times of the solution z(t,j) is locally finite. 

Proof: First of all we remark that for each agent j, each time a switching occurs, there 
must elapse an interval of time of length at least — (|| g (o o)||. +A/2) before the agent i 
switches again: that is, for each agent the inter-switching intervals have lengths which are 
bounded from below. Indeed, the agent's state Xi(t,j) evolves with speed not larger than 
11-^1 loo (I l<7(0; 0)| |oo + A/2) (see Lemma 6) and, after each switching, it has to cover a distance 
of length at least A/2 before it fulfills a switching condition again. 

Let us now suppose that there exists a solution z(t,j) and a sequence of switching times 
tj, with (tj,j) £ domz, such that j — > oo and 

j 

lim V(t fc+1 - t fc ) = T < +oo . (15) 

Because the number of switches is infinite and the number of agents is finite, then there is 
at least an index i £ V for which switches infinitely many times. We can then extract the 
sequence U k of times at which qi changes its value. Clearly, 

j 

lim V(t lfc+1 -t ik ) < +oo . 

i->+oo * — ' 
fc=0 

On the other hand, thanks to the previous argument, one has that 

A/2 

Uk+1 _ Uk ~ l|£||oo(||*7(0,0)|U + A/2)' 



and then 



3 

lim V(t 4fc+1 -t ik ) 

k=0 



i.e. a contradiction. 



Remark 3 (On the "continuous component" of the hybrid time domain) Combining 
Lemma 5 and Proposition 5 it is possible to conclude that the continuous component of the 
hybrid time domain is unbounded. Namely domz is either the union of infinitely many 
intervals of the form \tj,tj + {\ x {j}, with ~ to < ti < . . ., or the union of finitely many of 
such intervals, the last one being [tj, +oo) U {j}. Moreover, in the former case, necessarily 
Uj^[tj,tj + i] = R>o by Proposition 5. 
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Remark 4 (Data rate) In the case of quantizers with hysteresis, thanks to Lemma 6 and 
Proposition 5, one can find an upper bound on the the data rate at which each agent 
transmits information to its neighbors. In fact, the proof of Proposition 5 tells us that the 
information regarding the state is transmitted at most every T 



A/2 



l|i||oo(||g(0,0)|U+A/2) units 

of time. The information consists of packets of bits which encode the quantized state of the 
agents. In view of Lemma 6, the number of quantization levels employed by each agent is 



l|g(Q,Q)lh 



= +A / 2 +1. Hence, to encode these quantization levels, B := [log 2 ^8 



bits are needed. We conclude that each agent transmits information at a data rate which is 
not larger than 



B 
T 



log2(8 ^u +5 



J|g(0,0)lloo , 1 



\L\ 



We can also prove that the solution to (14) preserves the average of the continuous states. 

Lemma 7 (Average preservation) For each z(0,0) £ Xq, the solution z(t,j) to (14) is 
such that N- 1 l*x{t,j) = iV" 1 !*^;^, 0), for all (t,j) £ domz. 

Proof: By (13), for all (t,j) £ domz such that (t,j + 1) £ domz, x(t,j) = x(t,j + 1), 
i.e. during discrete transitions the value of x does not change and therefore the average is 
trivially preserved. For z(t,j) £ C, the continuous state x(t,j) satisfies N~ 1 l*x(t, j) — 
—N~ 1 l*Lq(t,j) for all t £ Ij. Since the graph is weight-balanced, then 1*L = 0*, and 
therefore the average is preserved also during continuous flow. This leads to the thesis. ■ 



Consider now the set of equilibria of the system, i.e. the set of states such that all the 
evolutions originating from those remain in the same state. This set is characterized in the 
following statement. 

Lemma 8 (Equilibria) The set of equilibria for the system (14) is £ = £\ U £2, with 

£1 = {(x, q) £ X : 3k £ Z s.t. q t = kA, and - ^ + < Xi < - + q i} Mi £ V} 

£2 = {(x, q) £ X : 3k £ Z s.t. q t = kA + § , and -| + % <Xj<f + qi, Vi £ V} . 

Proof: Suppose that the initial condition z(0,0) belongs to £\. Then, since g,(0,0) = 
fcA and -f + %(0,0) < 2^(0,0) < f + ft(0,0) for all i £ V, z(0,0) £ C, no discrete 
transition is triggered, and moreover —Lq(0, 0) = 0, i.e. during the continuous flow the state 
remains unchanged. The same can be argued if the initial condition z(0,0) belongs to £2- 

Suppose now that z(0, 0) £ X is an equilibrium point, i.e. z(t,j) — z(0, 0) for all (t,j) £ 
domz. Then z(0, 0) ^ D, because otherwise a discrete update via (13) would occur leading 
to q(0, 1) q(0, 0), a contradiction. Since X = C U D, then z(0, 0) £ C. In order to have 
x(t, 0) = x(0,0) for all t £ Iq, it must necessarily be true that — Lq(0,0) = 0, i.e. g(0,0) £ 
spanin(AZU(AZ + A/2)) JV or, equivalently, that for each i £ V, either ^(0,0) = fcA + f or 
9i(0,0) = kA for some k £ Z. Since z(0, 0) £ C, then necessarily +<&(0, 0) < x,;(0,0) < 
f + qi(0, 0), for all i £ V. We conclude that z(0, 0) £ £ x U £ 2 - ■ 

Differently from the case of quantizers with no hysteresis, the set of equilibria £ may not 
be globally attractive for the solutions to (14) as the following example shows. 
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(16) 



Example 2 (Finite-time limit cycle) Consider the system (12)-(13) with N = 2 and 
where x = —Lq is 

X\ = <?2 - qi 
Let qi G AZ be fixed and consider the initial condition 

zi(0, 0) = ( qi - p gi), 22(0, 0) = {qi + ^,qi + A). 



The continuous dynamics is given by 



±i = A 
±2 = -A 



so that 



xi(t) = gi-4 + At 



x 2 (t) = qi + fy-At 



on the interval [0, j] x {0}. Indeed, at time t = \ we have Xi(^) = gi, 2:2(3) — 9i + y> so 
that a discrete transition is triggered. It holds 

zi^' 1 ) = Z2 (^' 1 ) = (Si + + 2")- 

The continuous dynamics is then given by 

±1 = f 
i 2 = -f 

so that 

Xl(t) = g 1 + f(t-I) 

a*(t) = 9i + #-#(*- |) 

on the interval [j, |] x {1}. At time t = | we have xi(|) = 91 + -yi x 2(§) = 91, so that a 
new discrete transition is triggered. We have: 

.5 A A, ,5 „, . 

zi(-,2) = ( qi + y,9i + J"), ^(4,2) = (91, <7i)- 

In the next interval the continuous dynamics is given by 

*i = -# 
±2 = # 

so that 

11 (t) = 9i + f-f(t-|) 

i 2 (t) = 31 + #(*-!)■ 

At time t = | we have xi(|) = = £1(2) and x 2 (j) = qi + — = £2(3), so that a new 
discrete transition is triggered. Hence, z(|,3) = z(4, 1) and the evolution has entered a 
cycle. 

The conclusion is that there exist trajectories of system (12)-(13) which converge to 
periodic trajectories in finite time, and in particular do not converge to the set of equilibria. 
In this bidimensional example, the exhibited periodic trajectory lies in the closure of the set 
of equilibria. I 
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It should be noted that finite-time convergence of trajectories implies that solutions to (14) 
are not backward unique, whereas forward uniqueness has been proved in Lemma 5. 

In view of the existence of limit cycles, we can prove a convergence result analogous to 
Corollary 1. 

Theorem 3 (Finite-time convergence) Consider the system (12)-(13). Then, for any 
e G (0, 1), there exists T(e) > such that the solution z(t,j) to (12)-(13) satisfies 

Mtj)- ^,(0,0)11 <^ e -^M^*VN 

for all (t,j) G domz such that t > T(e). 

Proof: The proof follows the lines of those for Theorem 1 and Corollary 1. Define 
n* 

N 



y(t,j) = Qx(t,j) = (I — ^h-)x(t,j) and observe that for z(t,j) G C, y(t,j) satisfies 



y(t,j) = (I-^)x(t,j) = -(I-^)Lq(t,j) = -Lq(t,j), 

while at each switching time U, y(ti,i + 1) = y(ti,i). Given V(y) = we investigate 
dVlyV dt ^ = when z(t,j) G C. By Lemma 1, we obtain: 

VV{y){-Lq) = VV(y)(-L(x + q - x)) 

= -\7V(y)Lx-VV(y)L(q-x) 
= -VV(y)Ly-VV(y)L{q-x) 

< -\ 2 (Sym(L))\\y\\2 + \\y\\\\L\\^VN 



where in the inequality we exploit the fact that z G C implies q t — ^ < x t < q^ + and 
A 

2 



hence \xt — q%\ < 4 for all i G V. If 



IMI> l^ A 2 (SymU ^' (17) 

then 

VV(y)(-Lq) < - £ A 2 (Sym(L))|M| 2 = -2s\ 2 (Sym(L))V(y) . 
Hence, for all y(t,j) G C which satisfy (17), it holds 

At each switching time (U,i), on the other hand, we have V(y(ti,i + 1)) = V{y{t il i)). As a 
consequence, denoted by t = 0, t\, t 2 , ■ ■ ■ , tj the switching times which precede t, we have 

V(y(t,j)) < e-^^MW-VVfrfaJ)) 

< e -2eX 2 (Sym(L))(t-tj) V ^ tj j _ ^ 

< e -2cA 2 (Sym(L))(t-^) e -2eA 2 (Sym(X,))(^-^_i)y^^ ^ j _ -g^ 

< c -2eA 2 (Sym(L))(t-t J ) e -2eA 2 (Sym(L))(t 3 -t J _ 1 ) _ _ _ e -2eA 2 (Sym(L))(ti-i )^( 2/ ^ 0j Q\\ 
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which implies 

V(y(t,j)) < e^^^^V^O)) 

and therefore 

From the latter we conclude that, if at time (io,0), condition (17) is satisfied, then after at 
most 



T(£) = ^ \ ' eX 2 (SyL(L)) ln j } (18) 

units of time, ||y(i,j)||, with t > T(e), satisfies 

1 \\L\\ A /— 

lM < l-eA 2 (Sym(L)) 2 ' 

and continue to do so from that time on. The thesis then follows, recalling the definition of 
y(t,j) and that thanks to Lemma 7, l*x(t,j) = l*x(0,0) for all (t,j). ■ 



One might wonder whether it is possible to make a stronger claim of convergence to a 
set which does not depend on the network topology, for instance to the closure of the set of 
equilibria. The answer seems to be negative in general, as we show in the next subsection 
for systems with more than two agents by means of simulation examples: limit cycles are 
inherent to the hysteretic dynamics. An interesting open question, which we touch upon 
commenting the simulations, is to relate these limit cycles, and in particular their size, to 
the graph topology. 

Remark 5 (Convergence rate) From the proofs of Corollary 1 and of Theorem 3 we can 
see that, away from a strip containing the equilibria, the rate of convergence can be bounded 
by the same quantity for both algorithms. Actually, the rate is proportional to A2(Sym(L)), 
which is known to be the rate of convergence of the non-quantized linear consensus dynamics. 
In other words, our results suggest that the use of quantized measurements, with or without 
hysteresis, does not affect the rate of convergence of the consensus dynamics. 



4.3 Simulations 

In this section, we collect a few simulation examples, focusing on the effects of hysteresis and 
on the convergence properties of (14). We performed our simulations using a Matlab pro- 
gram which implements the hybrid system (14), and solves systems (5) and (12) by an explicit 
Euler scheme with time step St = 0.005. In the set of simulations we show here, we assume 
that N = 10, Aij € {0, 1} for every i, j £ V, and A = 0.05, and that the initial condition is 
x(0) = [0.91728, 0.26898, 0.76538, 0.18858, 0.28738, 0.09098, 0.57608, 0.68328, 0.54648, 0.42558]*. 
Let us first consider the two networks depicted in Figure 4. Figures 5 and 6 show for these 
networks the computed evolutions of (5) and (14), in terms of the quantized states q(x(t)) 
and q(t,j). Notice that solutions to (5) suffer from chattering. Indeed, there are subgraphs 
of the graph in Figure 4(b) which coincide with the graph discussed in Example 1. Hence, 
the chattering observed in the simulations is due to a sliding mode similar to the one de- 
scribed in the example and to the numerical implementation of the algorithm. Instead, 
solutions to (14) show no chattering: as expected, hysteresis prevents chattering, and may 
affect the convergence rate near the equilibria. The latter statement is not in contrast with 
Remark 5. In both the examples above, the quantized states of system (14) converge to 
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(a) 



(b) 



Figure 4: Two sample graphs, obtained as realizations of a random geometric graph [21] in 
the unit square [0, l] 2 with connectivity radius equal to 0.2. 




Figure 5: Examples of evolutions of quantized dynamics on the graph in Figure 4(a). 



consensus. To provide an example of limit cycles, we consider a directed ring topology, 
that is a graph such that A/j = {i + 1} for i e {1, . . . , N — 1} and JVn = 1. Figure 7 shows 
that in this case hysteresis induces a limit cycle of amplitude 2A on the quantized states. 
Furthermore, simulations demonstrate that larger cycles can be obtained on larger rings. 



5 Conclusion 

In this paper we analyzed a consensus problem subject to quantized communication, in 
terms of the stabilization of a dynamical system by means of a discontinuous feedback. Two 
quantization rules have been considered: a uniform static quantizer, and a hysteretic quan- 
tizer, which we designed with the goal of avoiding chattering phenomena. The convergence 
properties of both the resulting systems have been studied, and illustrated through simula- 
tions. Future work should include the analysis of other feedback consensus dynamics with 
quantization effects, and the application of these techniques to the problem of guaranteeing 
stable flocking of autonomous vehicles via quantized feedback control. 
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No hysteresis - quantized states 




Hysteresis - quantized states 




Figure 6: Examples of evolutions of quantized dynamics on the graph in Figure 4(b). 



No hysteresis - quantized states 





Hysteresis - quantized states 



Figure 7: Examples of evolutions of quantized dynamics on a directed ring. 
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